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1 Introduction 

It is well known that if a black hole is initially perturbed, the surrounding geometry 
will start vibrating into quasi-normal oscillation modes, whose frequencies and decay times 
depend only on the intrinsic features of the black hole itself, being insensitive to the details 
of the initial perturbation. The radiation associated with these modes is expected to be 
seen with gravitational wave detectors in the coming years, giving valuable information on 
the properties of black holes. For these reasons, quasi-normal modes (QNMs) of black holes 
in asymptotically flat spacetimes have been extensively studied (for reviews, see [Il[2]). 

The Anti-de Sitter - conformal field theory (AdS/CFT) correspondence has led to an 
intensive investigation of black hole QNMs in asymptotically AdS spacetimes. Quasi- 
normal modes in AdS spacetime were first computed for a conformally invariant scalar 
field, whose asymptotic behaviour is similar to fiat spacetime |3J. Subsequently, motivated 
by the AdS/CFT correspondence, Horowitz and Hubeny made a systematic computation 
of QNMs for scalar perturbations of Schwarzschild-AdS (S-AdS) spacetimes [1]. Their 
work was extended to electromagnetic and gravitational perturbations of S-AdS black holes 
in [5|. The study of scalar perturbations was further extended to the case of Reissner- 
Nordstrom-AdS (RN-AdS) black holes in [6]. Finally, the QNMs of scalar, electromagnetic 
and gravitational perturbations of RN-AdS black holes were presented in fT] using the 
results of [8]. 

The QNMs of AdS black holes have an interpretation in terms of the dual conformal 
field theory (CFT) According to the AdS/CFT correspondence, a large static black 
hole in AdS corresponds to an (approximately) thermal state in the CFT. Perturbing the 
black hole corresponds to perturbing this thermal state, and the decay of the perturbation 
describes the return to thermal equilibrium. So we obtain a prediction for the thermaliza- 
tion timescale in the strongly coupled CFT. In ref. [1] it was shown that the QNMs for 
the scalar perturbations of large Schwarzschild-AdS black holes scaled with the tempera- 
ture and it was argued that the perturbed system in the dual description will approach to 
thermal equilibrium of the boundary conformal field theory. However, when the black hole 
size is comparable to the AdS length scale there is a clear departure from this behaviour. 
It was then conjectured that this behaviour may be connected with a Hawking-Page phase 
transition pTl [T2] which occurs when the temperature lowers. 

These results were further confirmed in El El E] • However, the behaviour of QNMs 
for small black holes is still poorly understood. Another interesting finding of the elec- 
tromagnetic and gravitational perturbations is that purely dissipative modes appear in 
the spectrum which are pure imaginary QNMs. In such perturbed classical backgrounds 
the presence of these dissipative modes indicate that the boundary theory reaches ther- 
mal equilibrium with no oscillations. It was shown in ^ that for axial perturbations of 
Schwarzschild-AdS black holes these highly damped modes scale as the inverse of the black 
hole radius and this behaviour persisted in the case of Reissner- Nordstrom- AdS black holes 

It was observed in [T3] that for small black holes as we lower the temperature to a critical 
value there is a phase transition of a vacuum topological black hole towards a hairy black 
hole (MTZ). This claim was supported in [2] by calculating the QNMs of electromagnetic 
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perturbations of the MTZ black hole and topological black holes. It was found that there 
is a change in the slope of the QNMs as we decrease the value of the horizon radius below 
a critical value, and this change has been attributed to the phase transition. It was also 
observed in \1A\ that for small black holes the quasi-normal frequencies converge toward the 
imaginary axis, i.e., their real part decreases and after the first few quasi-normal frequencies, 
it vanishes, indicating that for black hole radius smaller than the AdS length scale there 
are only a finite number of QNMs. It was shown that the finite number of such modes 
for small horizons is due to the existence of bound states behind the horizon, which is an 
unobservable region. Further evidence that the behaviour of QNMs may provide indications 
of a phase transition was provided in |15j . 

Topological black holes [16]- [25], having hyperbolic horizons, introduce new features 
not present in spherical black holes. In [26] it was shown that hyperbolic black holes can 
be described as thermal Rindler states of the dual conformal field theory in flat space. It 
was also found that, for small topological black holes of the size of the AdS length scale, 
the entropy at strong coupling is larger than the entropy obtained from field theory at 
lowest perturbation order at weak coupling. One possible explanation put forward was 
that there is a phase transition of topological black hole to vacuum AdS space at a critical 
temperature, which however was not observed in |26j (see also |2Dl 127] ). 

In this work we make a detailed study of electromagnetic and gravitational perturbations 
of charged topological black holes (CTBH) in AdS space. Studying these perturbations of 
the background geometry, we show that the second-order phase transition observed at a 
critical temperature in [13] occurs in more general configurations including charge. We 
calculate both analytically and numerically the QNMs of axial and in some cases the polar 
perturbations of large and small black holes. We find that for large black holes the QNMs 
irrespectively of the value of the charge, exhibit a negative slope. However, if the value 
of the black hole radius is smaller than the AdS length scale and the charge is small, the 
propagating QNMs (whose frequencies have a non-vanishing real part) of both axial and 
polar perturbations are finite with a positive slope. We attribute this behaviour to a second 
order phase transition which occurs as the temperature approaches a critical value (or the 
black hole radius is approaching the length scale of the AdS space). We also find that, in 
the case of small black holes, if we increase the charge, the number of propagating QNMs 
is again infinite and a part of positive slope coexists with the negative slope frequencies. 
This observation indicates that in the case of small black holes, the charge plays the role 
of an order parameter and according to the AdS/CFT correspondence we expect that the 
dual boundary theory is described by a thermal state with coexisting phases. 

We also study the purely dissipative modes appearing in the spectrum, both analyt- 
ically and numerically of electromagnetic and gravitational perturbations of the charged 
topological black holes. For large black holes the purely dissipative modes of both ax- 
ial and polar perturbations scale linearly with temperature. Also the intermediate black 
holes to a good approximation depend linearly on the temperature. This is the expected 
behaviour for large and intermediate black holes and they agree with the normal QNMs 
behaviour discussed in [3]. If the black hole radius is smaller than the AdS length scale 
then we find a clear departure from linearity with temperature. If the charge is small the 
purely dissipative modes scale with the temperature according to a + 6/(T — Tq) where a, b 
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are constants. Then, for a fixed cliarge to mass ratio we observed an infinite change of 
slope at T = To signaling a second order phase transition. As we increase the charge the 
temperature dependence changes drastically. 

The paper is organized as follows. In section [2] we review the basic properties of the 
topological black holes and we discuss their thermodynamics. In section [3] we present the 
analytical calculations of the QNMs of electromagnetic and gravitational perturbations of 
the CTBH. In section H] we study numerically their behaviour and in section [5] we discuss in 
detail the purely dissipative modes appearing in the spectrum. Finally, section [HI contains 
our summary. 



2 Thermodynamics 

We consider the action 



a X\/—g 



where I is the AdS radius. The presence of a negative cosmological constant (A = — 3/^^) 
allows the existence of black holes with a topology i?^ x S, where S is a two-dimensional 
manifold of constant curvature. These black holes are known as topological black holes 
[Tn]-I2n]- The simplest solution of this kind, when S has negative constant curvature, reads 

ds' = - f-1 + ^ - ^1 de + ,au\ + ' (2-1) 

where da^ is the line element of S, which is locally isomorphic to the hyperbolic manifold 
and S must be of the form 

S = H'^/T with r C 0(2, 1) , 

where F is a freely acting discrete subgroup (i.e., without fixed points). 

The configurations (12.11) are asymptotically locally AdS spacetimes. It has been shown 
in [21] that the massless configurations where S has negative constant curvature are stable 
under gravitational perturbations. More recently the stability of the topological black holes 
was discussed in [22] • 

If we introduce to the above action an electromagnetic field 



1- j d'x^gF^^^F,, 



IGtt 

then the metric of a charged topological black hole is given by 



ds' = -f{r)de + ^^+r'da\ /(r) = - 1 - ^ + ^ , (2.2) 
/(r) r 



where we have fixed the length of the AdS space to / = 1. The horizon r+ is specified from 

= 2Gfir+ + rl-ri . (2.3) 
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Define the charge to mass ratio by 



(2.4) 



Then usine; Q we hav(0 



l + ^/l-A2(r2 -1) 



(2.5) 



We consider an electric potential at the horizon $ = —qjr^ and the electric charge Q 
aq/AiT. The temperature, the entropy and the mass of the CTBH are given by 



T 

S 
M 



f'{r+) _ 3rX -rl-q^ 



2 



AG ' 

cr/i (y{j'X — r\ + q^) 
47r SvrGri 



With the help of 



dM 



^(3ri - 1 - ^) 
SttG 



dr+ + ^ ^ dq 



47rr I 



it is easy to verify that the law of thermodynamics 

dM = TdS - ^dQ 

is valid. The free energy is 



F = M -TS + 



(j(r^ + r+) crg^ 



levrG 167rr+ IGnG 
where Gfi should be replaced by its value in eq. (12.51) and then 

dF = -SdT + Qd^ 
is easily verified. Notice now that the critical temperature 

T - ^ 



^ There is another possibility, 



X^jGfif 
r+ 



(2.6) 



(2.7) 



(2i 



(2.9) 



(2.10) 



l-^l-A2(r^-l) 

We shall not discuss it here because it corresponds to configurations which do not approach the critical 
point r+ = 1. 
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corresponds to r+ = 1. We have already given the functions F{r^) and T(r+), so we can 
readily Taylor expand F{T) about Tq and express the necessary derivatives through the 
formulae 

dF _ F' (fF _ T'F" - F'T" d^F _ ^F'T"^ + T'^F^^) - ^.T'T'F" - F'T'T^^^ 

dT ~ T'dT^ ~ ' rfT3 ~ ' ' ' ' ' 

where primes denote differentiations with respect to r+. The expansion of the free energy 
finally reads: 



(l + 2(T - To)7r + (2 + X^){T - To) 



P±^^!i(i±^(r-T„)VH-...). (2.11) 



It is interesting to consider how formulae (12.31) - (12. lip simplify in two important special 
cases, namely when takes on the values —1 or 0. Notice that having chosen a definite 
value for A means that q and r_|_ are no longer independent, 
(a) A^ = -1 (CTBH with scalar hair [28j) 

ds^ = -f{r)df + ^ + r^da^ , /(r) =r'-(l + ^)\ (2.12) 



/(r) V r J 

The horizon r+ is specified from 



G'/. = r+(r+-l) , (2.13) 

and the relevant quantities become 

O = 

in 

f(r+) 2r+ 



Q^?? (2.14) 



T 



An 2tt 
ar'i_ . „ ar.dr. 



' - 2G ' 

M = (2.15) 

dM = TdS - ^dQ , (2.16) 

F = M-T5 + g$ = --^ , (2.17) 

dF = -SdT + Qd^ , (2.18) 



O" 



1 + 2(T - To)7r + (T - Tq) V + 0(T - Tc 



SttG" 

(b) A^ = (Uncharged topological black hole) with 



/(r) = r'-l-^, (2.19) 
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2G/i = - r+ , (2.20) 
Q = aq/A-K = . (2.21) 



T 



f{r^) _ 3rl-l 
4:71 47rr I 



46" 2G ' 

M = ^ = ^ = ^M_I) rf,^ . (2.22) 

dM = TdS , (2.23) 

F = M-T5 = -^^^i^^ , (2.24) 

(iF = -SdT , (2.25) 



cr 



1 + 2(r - ro)7r + 2(T - Tof7i^ + (T - To)^^ + ...). (2.26) 



SttG" 

As an application of the above considerations if we assume that 

2 



/(r) = - (l + ^) (2.27) 

solves a four dimensional gravitational action describing a gravitational field coupled to a 
scalar and to an EM field (charged MTZ black hole [28]), calculating the quantity 



7rV(2 + 3A^)( 1 + 
16G 



^CTBH - FcMTZ = 7 ' " ^ {T - To)V + ... (2.28) 



we see that it is changing sign as we cross the critical temperature, hence we have a phase 
transition of the charged MTZ black hole to CTBH. 

3 Gravitational and Electromagnetic Modes 

The radial wave equation for gravitational perturbations in the black-hole background 
can be cast into a Schrodinger-like form [29] . 

-^ + V[rin)]^ = 0^^ , (3.1) 

in terms of the tortoise coordinate defined by 

±1 = J__ (3.2) 

dr f{r) 
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The potential V is determined by the type of perturbation. For axial perturbations, we 
y(r) = (r) = fir) \ + — f ' (3-3) 



where 



A = Jl + ^(A + 2)A2 , (3.4) 



and A are the eigenvalues resulting from the the solution of Laplace-Beltrami equation 
Aijj = —Alp on a space which has constant negative curvature and which is of the form 
P^i_i_.^(coshr)e*''"'^(coshr) and the corresponding eigenvalues read: A = — /(/ + 1) with 
/ = -i ± i^, so that A = 1 + 

For polar perturbations we have, 



where 



?7±(r) = 18(1 ±A)2G'V + ^V4^(1±^)'^ + 9(1 ±^)'(5t3A)-^'^' 



2 A + 2 ' ' ' ^ ' ' 

^2, ,2 



+18(1 ± A)(A + 2)^-^ + 3(A + 2)2(1 ± A)^ + A(A + 2f . (3.6) 

We shall refer to the modes corresponding to the potential V~ (axial or polar) as Zi modes 
whereas those corresponding to will be referred to as Z2 modes. 
The potential for polar perturbations (13. 5p has a singularity at 

r = rn = : (3.7) 

" A+2 ^ ' 

in addition to the standard singularities r = r+, 00. Remarkably, at the singularity (13. 7p . 

f/±(ro) = 2(A + 2)V(ro) , (3.8) 
resulting in the behaviour of the wavefunction 

^ ~ (r-roT , a = 2,-1 . (3.9) 

It turns out that at the boundary (r 00), the wavefunction has the same behavior with 
a = — 1. Precisely, 

^~(l-ro/r)-^ (3.10) 

so it does not obey the Dirichlet boundary condition \1/ = 0. Instead, it obeys the Robin 
boundary condition [30l [3l] 

r2^' + ro^ = (r^cx)). (3.11) 
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This is true even in the limit g = for the mode. For the Zx mode, ro = if g = 0, 
so we obtain the standard \1/ — ^ boundary condition. The latter is the electromagnetic 
mode. Indeed, if g = 0, A = 1, so 



[/_(r) = (A + 2)2A , V_{r)=k 



fir) 



(3.12) 



Evidently, the potential vanishes at the horizon (K(r+) = 0, since f{r^) = 0). This is the 
case for all types of perturbation. 

To obtain analytic expressions for the quasi-normal frequencies, it is convenient to 
introduce the coordinate [31] 

u = —. (3.13) 
r 

The wave equation fl3.ll) becomes 



V 



^ = , iu 



where prime denotes differentiation with respect to u and we have defined 



(3.14) 



f{u) 



fir) 



M , Viu) = -— 

I M I 4- 



y2 J.2 ^ 



(3.15) 



For the various potentials, we obtain 



-U H 



k = — 



(3.16) 



and 



V^{u) 



u'f{u)U±{u) 



(A + 3(l± A)^m)2 



(3.17) 



UM = 



18(1 ± A)' 



,GV , 811-A%^ , .^2GV^4^9(l^^)2(5^3^)_ 



r 



^ ^ 2 A 



:i±A)^ 



-u 



r 



+ 



-18(1 ± A)A^^z2 ^ 3^2(^ ± A)S(^n + A^ - 2K^\ , 



(3.18) 



where 



A = A + ^ 



(3.19) 
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3.1 QNMs of Large Black Holes 

To study the form of quasi-normal modes for large horizons, it is convenient to factor out 
the behaviour of the wavefunction at the horizon {u = 1), 

= {l-uy'^F{u) . (3.20) 

The wave equation becomes 

Ar^F" + Br^^^F' + C^ f^F = , (3.21) 

where 



^t^Th I — u 



"•^ f f (47rTH)2 (1 - uy ^ AttTh l~u^ AnTn {1 - uy ^ ' ' 

For V = (axial Z2 modes), this wave equation may be solved for small ci), A. To employ 
perturbation theory, write eq. 03.211) as 

{no + ni)F = 0, (3.23) 

where 

HoF = AooF" + BoofiF' + CofiF , 

HiF = (Ar^-Aoo)F" + iBr^,^-B^,o)F' + iC^ i-Co,o)F . (3.24) 

The zeroth order equation HqFq = is obtained by letting ci) , A — 0, r+ 00 while 
keeping Gfi/r"^ and g^/r^ fixed. We have A — >■ 1, so 



iu'Mn)F^y 



-u H — 



Fo = . (3.25) 



Despite its apparent complexity, the acceptable solution takes a remarkably simple form, 

Fo{u) = u . (3.26) 

The first-order constraint reads 

FoHiFo = , (3.27) 

Jo 

which imposes a constraint on the parameters (dispersion relation) of the form 

ao — isLiu — a2Lj'^ = . (3.28) 
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After some algebra, we arrive at explicit expressions for the coefficients, 



ai 



- / 

10 



J 2 SGufA- 1) 
-A - ^ + — ^ -u 

r± rl 



^ + ^ ^ 3G/i(l- A) 



4:7rTH Jo 



u 



2ufo ^( u^fo 



1-u 



1-u 



u'fo 



AtxTh \l-u 



(3.29) 



The other coefficient, a2 is not needed for the lowest mode. Therefore, 

2 



lu = ao 



^ + ^ 3G/i(l - A) 



+ 



4^3 



Explicitly, for large r+. 



(A + 2)rl 
QGfj, 



(3.30) 



(3.31) 



which is a purely dissipative mode. Notice that u oc l/r+, because Gfi oc for large r+. 

For Zi axial perturbations {V = V^), there is no solution in the small a) limit. This 
indicates that the lowest lying modes are proportional to r_|_ (so that a) = w/r^ remains 
finite as r+ oo). Explicit analytic expressions cannot be obtained in general. However, 
as we discuss later, in the case of g = the wavefunction may be written in terms of a 
Heun function leading to semi-analytic expressions for the frequencies. 

The calculation of polar modes is considerably more involved due to the additional 
singularity of the potential [31j . The spectrum is similar to the spectrum of axial modes as 
evidenced by our numerical calculations (section H]). We shall not perform the analytical 
calculation of these modes for large black holes as our main focus is on the critical point 
(r+ = 1) to which we turn next. 



3.2 QNMs at the Critical Point 

At the critical point (r_|_ = 1, g = 0, = 0), the wave equation for all types of perturbations 
reduces to 



UJ 



A 



^ = 



(3.32) 



whose solutions can be written in terms of associated Legendre functions. The solution 
which is well-behaved at the horizon is 



^(u) = PI 



u 



To see that it is regular at m = 1, express it in terms of a hypergeometric function. 



/I 4-U\ 11 1 

L{l — tuj)\l — u/ 2 2 2 



-)• 



At the boundary, m — > 0, it approaches a constant. 



^(0)=i^i/2(0) 



2'-v^ 



r(! - l^e - h^Mi + K 



(3.33) 



(3.34) 



(3.35) 
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Demanding that it vanish, we deduce the quasi-normal frequencies 

3^ 



LJn = ±^-i(2n + -^ , n = 0,1,2,... 



(3.36) 



which have finite real part (except in the special case ^ = 0). 
For explicit expressions, use 



F(a, f3; 7; z) = {1 - z^-'^-^Fi^ _ «, 7 - /3; 7; z) 



to write 



These hypergeometric functions are polynomials. Explicitly, 



(3.37) 



2 + 2ze 2 
1 H — u 



(3.38) 



etc. They are orthogonal under the inner product (no complex conjugation!) 

"1 du 



[n\m) = 



1 — w 



(3.39) 



defined by appropriate analytic continuation of the parameter ^. To normalize them 
{{n\n) — 1), choose 



4r(-iO 



v^r(-ie - i) 



42 - 
) ^1 — 



6r(-ie - 1) 



- |)r(-ie - i) ' 



(3.40) 



etc. 



1 


{n\'H'\n) 


2uJn 


(n n) 



Moving away from the critical point, the frequencies shift by 

5uJn 

where 

where we applied standard first-order perturbation theory. 
We obtain for the axial modes 



(3.41) 
(3.42) 



6uJi 



i 1- 



+ + 



'3(1 ± A) 



4r(-20 
v^r(| - ^0 



1±A 



7 _ {f (l±A)-2 + 9(l + 
2 + 0Fr(| - 



e-i4e2}r(-i-ze) 



(3.43) 
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For small ^, the change in the imaginary part is negligible whereas the change in the real 
part is 

1 \ 3(1 ± A) 

For Zi modes, above the critical point (r+ > 1), 5a;o < and the real part decreases. There 
is a critical value of ^ (determined by ^ + boj^ ^ 0) , 



Sujo^-{l- — ] ' . (3.44) 



Co - V12(l + A)(T-ro) , (3.45) 

where we used T — Tq i^^l"^^' below which the mode does not propagate (purely 

dissipative mode). It turns out that for ^ < there is a pair of purely dissipative modes. 
The first harmonic behaves similarly with a higher critical value of ^, 



6~ V[27(l + A)-8](T-To) (3.46) 

below which it turns into a pair of purely dissipative modes. 

Below the critical point (r+ < 1), 6uJo > and the real part of the Z2 modes increases. 
The modes do not become purely dissipative for any value of ^. 

Also notice that above the critical point, ScUn increases with n, therefore the real part 
decreases with n (positive slope) whereas below the critical point we obtain a negative slope 
for propagating modes. 

Zi modes behave in the opposite way because 1 — A < 0. Above the critical point, 
these modes never become purely dissipative. Below the critical point, we obtain the 
critical values 



eo « V12(A - 1)(T - To) , 6 ~ + 27(A - 1)](T - To) (3.47) 

similar to Z2 modes above the critical point. 
For the polar modes we obtain 



Scuq — i (1 — 



6lji = ill \ 



7 . {-1(1 ± (1 + - 2 + 9^e - Ue] r(-l - 

2^'^ Ar(i - iO 



(3.48) 

These modes coincide with their axial counterparts for A = —1 (electromagnetic modes). 
3.3 Uncharged Black Holes 

In the case of no charge in the black hole, the wave equations simplify because v?f{u) has 
at most three distinct roots, u = 1,1], —j^, where 

V = ; ^ . (3.49) 

1 + ^1 - 4(1 - 1/rl) 
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Then the solution to the wave equation may be written in terms of a Heun function. 

For Zx perturbations, the potentials for axial and polar modes coincide reducing to the 
electromagnetic potential. The wavefunction may be written as 



1 —u 

X Heun(a, g, a, /3, 7, 5, ) . (3.50) 

1 — T] 

The Heun function obeys the equation 

z{z-l){z-a)F"-[{-a-i3-l)z^ + {{d+-f)a-d+a+p+l)z-a-f]F'-[-ai3z+q]F = (3.51) 
and the various constants are 



l + 2r) .1+77 + 77 



2 



q= -A 



1 oil -1 9 ' 

1—77^ 1 — 77"^ 

^ _ l + rj-2{l-too)7^^ 2 - {2tu; + l)r/ - 77^ 

' ^ ' ^ (l-r/)(l + 27/) ' (l-r;)(2 + 7;) " 

It behaves nicely at the horizon (li — > 1). Requiring \E'(0) = yields the constraint 

}ieun{a,q,a, P,j,S,- )=0, (3.52) 

which may be solved for uj to obtain the quasi-normal frequencies of axial Zi modes. 
As r_|_ — > 00, the two lowest purely dissipative modes asymptote respectively to 

3 

cuq — —i-r^ , Ui — — 3ir+ (3.53) 

if A is kept constant. As A increases, the two modes approach each other coalescing at 

A = 0.115r^ , uo = oui = -2mir+ (3.54) 

Beyond this point, they develop a finite real part. 

For axial Z2 perturbations, the solution to the wave equation may be similarly written 

as 



■ - J, ( ri \ ""(2+'7)(l+27j) ^2 

^(li) = (it - r^) (1-^X2+,,) U + (1 _ ii)«^(i-^)(i+2^) 

1 — n 

X Heun(a, g'', a', 7, 5, ), (3.55) 

1 — 7/ 

where the parameters a, 7, 5 are same as before and 

1 — 77"* 
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It behaves nicely at the horizon {u — > 1). Requiring \I/(0) = yields the constraint 

1 



Heun(a, q\ a', (3', 7, 6, 



1 — rj' 







(3.57) 



which may be solved for ci) to obtain the quasi-normal frequencies of axial Z2 modes. 
As — > CX3, we find a single purely dissipative axial Z2 mode which asymptotes to 



A + 2 



UJo = -i- 



3r I 



(3.58) 



confirming the earlier analytical result (13.311) in the case g = 0. 

The above results do not apply to polar Z2 perturbations due to the additional singu- 
larity in the potential which survives the limit g — > 0. 

4 Numerical Calculations and Results 

We briefly review the method of Horowitz and Hubeny [1] as it is applied to our problem. 
After performing the transformation \l'(r) = il)^{r)e~'^'^^* , the wave equation (13. ip becomes 



(Pi)u{r) ( df{r) \ di)^{r) 



dr"^ 



dr 



dr 



(4.1) 



where the potential V{r) is given by (I3.16P or (I3.17p . The change of variables r = 1/x 
yields an equation of the form 



(X — Xj 



, d^4j^{x) 
dx'^ 



+ t(x) 



(X — X. 



dil)^{x) 
dx 



+ u{x)i)^{x) = 



where x+ = l/r+ and s{x),t{x) and u{x) are given by 

k 

^i^) = '^Sn{x-X+)'' 
n=0 
k 

n=0 
k 



uix] 



n=0 



where k = 3 for axial and k = 7 for polar perturbations. Expanding the wavefunction 
around the (inverse) horizon x^, 



i^U^) = E «n(t^)(a; - a;+)" 



(4.2) 



we arrive at a recurrence formula for the coefficients, 



a„(u;) 



1 



n-l 



n(n - l)so + nto + Uo 



J2 ["^("^ - 1)-S„,_m + mtn-m + M„-m]am(w) . (4.3) 



m=n— 3 
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Figure 1: QNMs of axial perturbations with r+ = 1.50, v = 0.10. 

We note that the few coefficients am{^^) with negative index m which will appear for n < 2 
should be set to zero, while ao(a;) is set to one. Since the wave function should vanish at 
infinity (r — oo, x = 0), we deduce 

oo 

V;^(0)=^a„H(-x+r = 0. (4.4) 



The solutions of this equation are precisely the quasi-normal frequencies. 
4.1 Overview and Spacing 

3^4 _q2 

To begin with, let us recall the expression for the temperature given earlier: T = """^^^t — . 

The constraint that T should be positive yields the inequality < q^^x = — r'^, so the 
charge may be expressed through the parameter 



z/^-^ , g^.. = f v'3<-^+- (4-5) 

We note that this parameter is related to the charge-to-mass ratio A through the equation 

4z/2(3r2 - 1) 



A^ 



{3u'^r'^ + r\ — u"^ — 1)^ 



One needs to examine separately large (r+ > 1) and small (r+ < 1) horizons. A typical 
graph for axial perturbations at r+ = 1.50, z/ = 0.10 (A = 0.37) may be seen in Fig. [TJ 
The figure depicts the curves 3fJ['?/'(ci;)] = and ^^[^/'(ti;)] = in the complex u plane. The 
QNMs are given by the intersections of the curves. We remark that, if one views uji versus 
ojR, the slope is negative, as can be seen in Fig. [U this feature does not change if one 
varies u. This specific graph depicts axial Zi QNMs at r+ = 1.50 and v = 0.10, but no 
qualitative change occurs if one considers Z2 rather than Zi perturbations. The behaviour 
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Figure 2: QNMs of axial perturbations with r+ = 0.95, v = 0.10. 

is similar for polar perturbations: for this value of r_|_ only some quantitative changes have 
been observed between the polar perturbations and their axial counterparts. 

However, for the small horizon r+ = 0.95 the behaviour of QNMs is different. For u = 
0.10 (A = 3.25) and axial Z2 perturbations (Fig. [2]) we find a finite number of propagating 
QNMs with a positive slope. (We note that there might also exist purely dissipative modes 
with vanishing real part, which are not clearly visible on such graphs; we will examine them 
in the next section). 

For u = 0.40 (A = 5.95) the pattern is qualitatively different as can be seen in Fig. [31 
a part with positive slope coexists with the negative slope frequencies and the number of 
QNMs is infinite again. We find similar results for the axial Zi perturbations, as well as 
for the polar Zi and Z2 perturbations. 

We may also make some quantitative statements about these modes. The asymptotic 
spacing between the QNMs may be easily read off from the figures and the results are 
presented in Table [H For = 0.95, the spacing is not constant; thus we report the 
difference between the lowest and the second lowest QNMs just to get some feeling for the 
order of magnitude. An important observation is that positive and negative slope parts 
coexist for r+ = 0.95 when u is large enough (Fig. [3]), so the results for the imaginary parts 
refer rather to the absolute value than to values with a well defined sign all the way. 

4.2 Lowest Modes 

For r+ > 1 the QNMs of the form presented in the previous section may be represented 
as a multiple of the spacing plus their lowest possible value, referred to as the offset in the 
literature. It is interesting to examine the behaviour of these lowest modes, in addition 
to the spacing. The charge is given as a fraction u of its maximal value and notice that 
we cannot increase the charge parameter u beyond a value about 0.40, due to convergence 
problems. 
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r+ 


V 


type 


^1 


Z2 


1.50 


0.10 


axial 


1.25 - 3.25i 


1.25 - 3.30i 


1.50 


0.30 


axial 


1.05 - 3.85i 


1.05 - 3.85i 


1.50 


0.00 


polar 


1.30- 3.40i 


1.36 - 3.30i 


1.50 


0.40 


polar 


0.94-4.20i 


1.03-4.17i 


0.95 


0.10 


axial 


0.26 + 1.67i 


0.09 + 2.01i 


0.95 


0.40 


axial 


0.00 + 1.80i 


0.17 + 2.20i 


0.95 


0.00 


polar 


0.32 + 1.80z 


0.38 + 1.70i 


0.95 


0.40 


polar 


1.90 + 2.30Z 


1.42 + 2.43i 



Table 1: Spacing of QNMs for various values of the parameters r+ and v with = 1. 

4.2.1 Axial Perturbations 

We examine the intermediate horizon r+ = 1.50 and a typical large horizon, namely r+ = 
20.00. It turns out that there are qualitative differences between the two. In Fig. [Hone may 
see the real parts of the axial Z\ QNMs for r+ = 20.00 for large horizons and small values 
of the charge purely dissipative modes are present, which is not the case for intermediate 
horizons, such as r+ = 1.50. We note that no different behaviour of this kind shows up for 
Z2 modes. 

In Fig. [5] we depict the absolute real part of the lowest axial modes for r_|_ = 0.95. We 
observe that, as one increases the charge, |3fJa;| for Z2 modes approaches a constant value, 
which is presumably also the value for the extremal black holes (z/ = 1); on the other hand, 
I^Rcul for Zx modes approaches zero. 

The agreement with the analytical result (13.431) is good (improving as we approach the 
critical point r+ — ^ 1). In Figs. [6] and [71 we show plots of the analytical expressions for 
|3fJu;| and Q'a;, respectively, for Z^ and Z2 axial perturbations for two different values of 
the horizon just below the critical point (r_|_ = 0.95 and r+ = 0.995). The accuracy of the 
analytical approximation decreases as the charge increases and fails when the first-order 
correction becomes comparable to the zeroth order approximation (I3.36p . Naturally, the 
point of failure increases as we approach the critical point. 

4.2.2 Polar Perturbations 

Qualitative differences show up for the small value r+ = 0.95, so we concentrate on this 
case. In Fig. [8] we present the real and imaginary parts of the axial and polar perturbations 
of Z\ and Z2. 

For Z2 polar perturbations the real parts (Fig. [HI left panel) do not have the striking 
behaviour of the axial perturbations, which are zero above v = 0.25. Similarly the imaginary 
parts (right panel) decrease with u with no spectacular change at z/ = 0.25, in contrast with 
their axial counterparts. Here we also observe the very mild behaviour of the polar QNMs 
with the charge as compared with the axial perturbations. 
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Figure 3: QNMs for axial Z2 perturbations with r+ = 0.95, u = 0.40. 

5 Purely Dissipative Modes 

After the overview of the propagating QNMs presented in the previous sections, we will 
make a systematic study of the purely dissipative modes of electromagnetic and gravita- 
tional perturbations. 

5.1 Axial Zi Modes 

Fig. [9] contains the results for purely dissipative QNMs at r_|_ = 5.00 and r_|_ = 1.50 versus u. 
For r_|_ = 5.00 and u = there are two purely dissipative modes with the smallest imaginary 
parts. For non-zero u these modes approach each other until finally, at z/ = 0.0663, they take 
on the same value and they acquire non- vanishing real parts, transforming into propagating 
QNMs. We obtain agreement with the asymptotic analytic expression f l3.53p for z/ = 
(^uj = -7.50, -15 for r+ = 5.00). 

For intermediate horizons, such as r_|_ = 1.50, the picture is different: the purely dissi- 
pative modes have relatively small absolute values for u = 0.40 and their absolute values 
increase as one approaches u = 0. We note that for u = there are no purely dissipative 
modes at all for this value of r+, so it appears that the imaginary parts of the QNMs tend to 
—00 as u ^ 0. This is in agreement with the modes obtained by solving the exact analytic 
equation (13.521) which is valid for z/ = 0. 

The QNMs for Zi perturbations fall into two disjoint classes; they belong to either large 
horizons (r+ > 3.15) or small ones (r+ < 3.15). They have different properties for = : 
large horizons have purely imaginary modes in this limit, while small ones do not. The 
behaviours of the two classes are qualitatively different. 
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Figure 4: The (absolute) real part of the axial Zi QNMs with r+ = 20.00 versus fractional 
charge u. The corresponding A ranges from to 0.05. The solid line corresponds to the real 
part of the lowest Zi mode; the dashed line is the real part of the second lowest Zi mode. 
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Figure 5: The (absolute) real parts of the axial Zi (lower curve) and Z2 lowest modes at 
r+ = 0.95 versus fractional charge u. 

If we go now to a horizon close to and below the critical point, r_|_ = 0.995 (Fig. [TOl) . we 
have in some sense the reverse behaviour: a finite number of propagating modes corresponds 
to the smallest absolute values of the imaginary part, a positive slope is observed and finally 
purely dissipative modes appear with large imaginary parts. For small u we observe the 
phenomenon just described, the propagating modes appearing last; for somewhat larger u 
purely dissipative and propagating QNMs are mixed. An important point is that the i/ — > 
limit of the imaginary parts is finite for r+ = 0.995, in contrast to what happens for the 
intermediate value r+ = 1.50 (Fig. M>- 

5.1.1 Temperature Dependence 

An issue that should be examined is the r+— dependence of the imaginary parts. It turns 
out that it is advantageous to use the temperature rather than r+ as a variable: in Fig. [11] 
we depict the (absolute value of the) imaginary part of the purely dissipative modes cor- 
responding to zero charge for large horizons versus the temperature. We observe that the 
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Figure 6: The (absolute) real parts of the axial Zi (lower curve) and Z2 lowest modes at 
r+ = 0.95 and r+ = 0.995 versus A using the analytical expressions (13.361) and (I3.43p . 




Figure 7: The imaginary parts of the axial Z\ (upper curve) and Z2 lowest modes at 
r_|_ = 0.95 and r+ = 0.995 versus A using the analytical expressions (13.361) and (I3.43p . 

imaginary parts scale linearly with the temperature to high accuracy. The data points have 
been calculated numerically as well as by solving the exact analytic equation (I3.52p . the 
two results being in excellent agreement with each other. The linear fits agree with the 
analytic asymptotic expression (13.531) . 

We also examined the intermediate horizons, in particular r+ = 0.90 up to = 2.00 
with A = 0.50. The relevant temperature T — Tq ranges from —0.04 to 0.24. The two lowest 
modes are depicted in Fig. [121 Good quality fits (xio./. < 1) of the form a + jSfr are 
possible for both branches: T > Tq and T < Tq. The curves for the lowest modes are also 
included in the figure. We observe that there is an (infinite!) change in the slope of this 
graph at T = Tq, which presumably signals a phase transition. 

In Fig. [131 we show a similar plot in the case of no charge but this time using the exact 
analytic eq. (13.531) . We obtain a similar singular behaviour {^uj ~ (T — Tq)~^) reinforcing 
the conclusion that a phase transition occurs at T = Tq. 
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Figure 8: Panel (a): The real parts of the lowest QNMs of the axial (lower curve) and polar 
(upper curve) perturbations for Z\ versus u at r+ = 0.95. Panel (b): The imaginary parts 
of the lowest QNMs of the axial (upper curve) and polar (lower curve) perturbations for 
Zi versus u at = 0.95. The real, panel (c) and imaginary, panel (d), parts of the lowest 
QNMs of the axial (upper curve) and polar (lower curve) perturbations for Z2 versus z/ at 
r+ = 0.95. 



5.1.2 ^—Dependence 

For large horizons (typical value r+ = 5.00) we consider the i/ = case. The ^—dependence 
is shown in Fig. [TH As ^ grows, the purely dissipative modes converge towards each other 
and finally disappear (that is, they turn into propagating modes with finite real part), a 
behaviour strongly reminiscent of the i/— dependence, depicted in Fig. [HI We have chosen to 
depict in Fig. [15] the ^—dependence of the lowest QNM for a typical intermediate horizon. 
Numerical results are in excellent agreement with the results obtained by solving the exact 
analytic equation fl3.52p . In particular, for ^ = 0, we obtain the two values '^u = —15, —7.50 
from eq. fl3.53l) whereas the two modes coalesce at the value = —10.25 at A = 2.875 
corresponding to ^ = 1.62 (eq. 03.541) ). all in agreement with numerical results. 

For r+ = 0.995, z/ = we know from the analytical calculations that the imaginary part 
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Figure 9: The imaginary part of the axial Z\ purely dissipative QNMs with r+ = 5.00 
(upper left corner), the corresponding A being smaller than 0.042 and r+ = 1.50 (right part 
of the figure) versus fractional charge the corresponding A ranging between 0.65 and 
0.90. 



of the QNMs changes little, so we focus on the ^—dependence of their real parts. We confirm 
numerically that the change in the imaginary parts is small and calculate numerically the 
real parts. The results for the lowest QNMs are shown in Fig. [ini Both numerical and 
analytical results (eq. fl3.43p ) are shown for comparison. 
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Figure 10: Im{u;) versus u for the five lowest Zi purely dissipative modes at r+ = 0.995. 
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Figure 11: The (absolute value of the) imaginary part of the two Zx purely dissipative 
QNMs with z/ = 0, ^ = 1 versus temperature. Horizons vary from r+ = 5.00 to r+ = 12.00. 
Also shown are the asymptotic expressions (13.531) . 



We observe that the real parts decrease with decreasing ^ until they vanish at some 
critical value. For ^ smaller than this critical value, the QNMs will be pure imaginary. 
The critical velues of ^ given by the analytic expressions (I3.47P are in good agreement with 
numerical results. This is seen in Fig. [17] which shows the singular behaviour of the three 
lowest critical values of ^ just below the critical point T = Tq. The data points have been 
calculated both numerically and by using the analytic equation (13.521) . 

For a given ^, the corresponding modes will be purely dissipative for large enough '^uj, 
since the relevant curve will have crossed the axis at a value greater than the given ^; 
however there will exist in general curves (that is, appropriate values of Qu), for which 
the corresponding ^uj will be different from zero and the QNM will be a propagating one. 
Therefore, QNMs with a sufficiently large absolute value of the imaginary part will be purely 
dissipative, while the lowest QNMs will be propagating. This calculation confirms once 
again that the lowest QNMs will be propagating, followed by purely dissipative ones. This 
fully agrees with the numerical results for small horizons and zero charge. The situation is 
exactly the opposite for large horizons: in that case the curves of approximately constant 
imaginary part never cross the horizontal axis. For non-zero values of the charge the picture 
is similar. 

5.2 Axial Z2 Perturbations 

The results are described starting in Fig. [TS] which provides a picture similar in some 
respects to the corresponding result for Zi excitations shown in Fig. O Two differences 
should be noted, however: (a) There is no qualitative distinction between intermediate and 
large horizons here. In particular, there is no need to also depict the behaviour of r+ = 5.00, 
as we did in Fig. [9], and (b) the absolute value of the lowest QNM is much smaller than 
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Figure 12: The imaginary parts of the two lowest axial purely dissipative QNMs with 
A = 0.50 versus temperature. For T < Tq the lowest modes are represented by points, 
while the second lowest ones with lines. For T > Tq the lowest modes are represented by 
points and lines, while the second lowest ones with points. Horizons vary from r_|_ = 0.90 
to r , = 2.00. 
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Figure 13: The imaginary part of the lowest frequency ujq of the axial Zi purely dissipative 
mode versus temperature for A = 0, ^ = 1 from the exact analytic eq. (I3.53p . Also shown 
is a fit Qujq = 
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Figure 14: The imaginary part of the two Z\ purely dissipative QNMs for r+ = 5.00, z/ = 
versus ^. 

the remaining ones; it almost coincides with the horizontal axis. 

For small horizons we have chosen to depict the results for r^. = 0.995 in Fig. [13 The 
picture is qualitatively the same as the corresponding result for Zx excitations, shown in 
Fig. [lOl In particular, nothing is special about the lowest mode in this case: this should 
be contrasted with the results above the critical point presented previously in Fig. [181 

5.2.1 Temperature Dependence 

Next we examine the temperature (or horizon) dependence of axial modes. Fig. [20] 
contains the numerical results for the two lowest purely dissipative modes versus T — Tq for 
A = 0.50. The lowest mode has been reported in [51 [7| to scale as ^ which agrees with our 
analytic asymptotic expression ( ]3.3ip . Using the temperature, rather than the horizon, as 

-14 I 1 1 1 , 1 
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Figure 15: The imaginary part of the axial lowest Zi purely dissipative QNMs with r+ = 
2.00 and u = 0.40 versus ^. 
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Figure 16: The real part of the four lowest axial Zi QNMs with r+ = 0.995 versus the ^ 
parameter for u = 0.00. The imaginary part decreases from left to right. The second panel 
shows the analytical estimate for the lowest two modes (eq. (13.43 1) ) for comparison. 

an independent variable, we have found that this mode scales as a+ jSf^ above the critical 
temperature; this fit is of very good quality and is shown in the figure. However, we have 
not been able to fit the T < Tq data to a function of this form except very close to the 
critical point T = Tq. There is an infinite change in slope at the critical temperature, as in 
the Zi case. We also remark that for T < Tq the two lowest values are very close to each 
other, while for T > Tq we have a lowest mode with a very small absolute value, while the 
absolute value of the second lowest mode is much larger. 

The above conclusions are confirmed by our analytic results. In Fig. [2T|, we show the 
lowest mode for A = for various temperatures both above and below the critical point. 
The data points were found by solving the analytic equation (13.571) . Once again, we observe 
a singular behaviour characterized by an infinite change in the slope at the critical point. 

5.2.2 ^—Dependence 

We now proceed with a discussion of the ^—dependence for the two regimes: large and 
small horizons. Fig. [22] contains both numerical and analytic results (using the asymptotic 
expression (13.311) ) for r+ = 20.00 The value ^ = yields the uppermost curve. This 
behaviour is similar to the behaviour for intermediate horizons in the axial Zi case (Fig. [T^ 
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Figure 17: The first three critical values of ^ for axial modes versus temperature for 
= 0. Also shown are fits, = 7r(To - T), = ./8(To-T), ^2 = a/20(To -T). 



and is quite different from the behaviour for large horizons in that case (Fig. [T4|) . 

For horizons below the critical point, results are contained in Fig. [23] where we chose 
a value close to the critical point, r+ = 0.995. The results are slightly different from the 
Zx case (Fig. [TB]) . In fact, the real part initially increases, it attains a maximal value and 
then it decreases and cuts the axis. This happens for very large values of the imaginary 
part, of the order of 30, the exact value depending on ^. Coming back to Fig. [23l the 
lowest plotted curve corresponds to such a value. We have also plotted the analytic results 
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Figure 18: versus fractional charge v for the five lowest Z^ purely dissipative modes at 
r , = 1.50. 
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Figure 19: The imaginary part of the axial purely dissipative QNMs at r+ = 0.995 
versus fractional charge v. 

f l3.36p and (13.431) . The agreement with the numerical results is good when ^ is not small. 
As ,^ — * 0, the corrections (13 .431) exceed the zeroth order result (I3.36p and the first-order 
approximation fails. 

For comparison, we also show results slightly above the critical point in Fig. [23] where 
we chose r+ = 1.005 and u = 0. Here the behaviour is very similar to the behaviour of Zi 
modes slightly below the critical points (c/. with Fig. [T6|l . Similar remarks can be made 
in this case. For each mode, there exists a critical value of ^ below which the mode does 
not propagate and becomes purely dissipative (Jituj = 0). These critical values of ^ depend 
on the temperature. Their behaviour near the critical tempetarure is given by the analytic 
expressions (13.451) . (13.461) and plotted in Fig. [25] (c/ with Fig. [T7]for Zi modes). Data points 
were found numerically as well as by solving the analytical equation (13.571) . 
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Figure 20: Quj versus T — Tq for the lowest and second lowest axial Z2 modes at A = 0.50. 
For T < Tq the lowest modes are represented by a line just connecting the points, while the 
second lowest ones with points. For T > Tq the lowest modes are represented by a fit and 
the points, while the second lowest ones with points. Horizon values vary from r+ = 0.90 
to 2.00. 
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Figure 21: The imaginary part of the lowest frequency ujq of the axial Z2 purely dissipative 
mode versus temperature for ^ = 1. Also shown is a fit '^Uq = ~^^:fZfo\- 

6 Conclusions 

We have studied the perturbative behaviour of the charged topological black holes. 
We have calculated both analytically and numerically the QNMs of electromagnetic and 
gravitational perturbations of these black holes. 

For large black holes we found analytically that the axial Z2 QNMs are purely dissi- 
pative, depend on the charge and scale as the inverse of the black hole horizon. The Zi 
axial modes are proportional to the black hole horizon but analytical expressions cannot 
be obtained in general for non-zero charge. For zero charge the potentials for both axial 
and polar modes reduce to electromagnetic potential and the wave function can be written 
in terms of the Heun function, leading to a semi-analytic expressions of the QNMs. 

For small black holes, at the critical point with zero charge and mass the wave equation 
simplifies and the QNMs can be explicitly calculated. For small changes around the critical 
point, the real part of the Z2 modes increases above the critical point giving a positive 
slope, whereas below the critical point it gives a negative slope. Above the critical point, 
for these modes there is a critical value of ^ below which there are only purely dissipative 
modes. Below the critical point, there are no purely dissipative modes for any value of ^. 
The Zi modes exhibit the opposite behaviour. Similar behaviour is exhibited by the polar 
modes. 

These results are also supported by numerical investigations of the QNMs. The numer- 
ical results show clearly a change of slope of QNMs around a critical temperature for all 
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Figure 22: '^uj versus fractional charge u for the lowest purely dissipative mode of axial 
Z2 perturbations at = 20.00. The corresponding A is smaller than 0.15. The curves 
correspond to (top to bottom) ^ = 0,1,2,3,4,5. The second panel shows the analytical 
expression for comparison. 

kinds of perturbations. We found that the purely dissipative modes scale linearly with tem- 
perature for large black holes, while for small horizons they scale according to a+6/(T--To). 
Then, for a fixed charge to mass ratio we observed an infinite change of slope at T = Tq 
signaling a second order phase transition. 

The numerical results show also an interesting dependence of the modes on the charge 
of the black hole. For small horizons and small charge the number of propagating {^u 7^ 
0) QNMs is finite, while as the charge increases, positive slope frequencies coexist with 
frequencies of negative slope and the number of propagating QNMs is again infinite. As 
the charge increases, a drastic change in the temperature dependence occurs. 



Figure 23: The (absolute values of the) real parts of modes versus ^ at r+ = 0.995. The 
curves correspond to the values (top to bottom) '^uj = —5.45, —3.47, —1.49, —29.21. The 
second panel shows the analytic approximation for the lowest two modes. 
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